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We numerically investigate dynamic critical behaviors of two-dimensional (2D) Josephson-j unction 
arrays with positional disorder in the scheme of the resistively shunted junction dynamics. Large- 
scale computation of the current voltage characteristics reveals an evidence supporting that a phase 
transition occurs at a nonzero critical temperature in the strong disorder regime, as well as in the 
weak disorder regime. The phase transition at weak disorder appears to belong to the Berezinskii- 
Kosterlitz-Thouless (BKT) type. In contrast, evidence for a non-BKT transition is found in the 
strong disorder regime. These results are consistent with the recent experiment on positionally 
disordered Josephson-j unction arrays; in particular, the critical temperature of the non-BKT transi- 
tion (ranging from 0.265 down to the minimum 0.22 in units of Ej/ks with the Josephson coupling 
strength Ej), the correlation length critical exponent v = 1.2, and the dynamic critical exponent 
z — 2.0 in the strong disorder regime agree with the existing studies of the 2D gauge-glass model. 



PACS numbers: 74.78.-w, 74.25.Fy, 74.81.Fa 
I. INTRODUCTION 

Phase transitions for the two-dimensional (2D) 
Josephson-junction arrays (JJAs) with weak positional 
disorder (i.e., on a slightly disordered lattice) and in a 
perpendicular magnetic field have attracted much atten- 
tion. Here positional disorder in the presence of an exter- 
nal magnetic field effectively induces random phase shifts 
of the magnetic bond anglesji*^ thus the corresponding 
2D random gauge XY (RGXY) model, in which the 
magnetic bond angles are quenched random variables dis- 
tributed in a certain width, provides a theoretical realiza- 
tion of such a positionally disordered Josephson-junction 
array (PDJJA) in a magnetic field: Strong magnetic 
field in the PDJJA corresponds to strong disorder in the 
RGXY model. It has been observed that sufficiently weak 
disorder in the KGXY model does not destroy the quasi- 
long-range order present at low temperatures and ac- 
cordingly, the system undergoes a Berezinskii-Kosterlitz- 
Thouless (BKT) type transition at a finite critical tem- 
perature Tc, which decreases as the disorder strength is 
raised up to a critical valuei^i^ On the other hand, in 
spite of a number of studies, the strong disorder regime 
of the RGXY model (and its fully disordered limit corre- 
sponding to the gauge-glass model) has resisted adequate 
understanding. In parallel with the Mermin- Wagner the- 
orem^ for the absence of the long-range order in the XY 
model, Nishimor:^ has proven the absence of the long- 
range glass order in the gauge-glass model. However, 
it should be noted that the vanishing local glass order 
parameter in the gauge glass is not completely incom- 



patible with the existence of the superconducting phase 
at nonzero temperatures. For instance, quasi-long-range 
glass order^ and possibility of a continuous phase transi- 
tion of an anomalous dimension^ have been suggested. 

On the one hand, there exist numerical evidences sup- 
porting the zero-temperature phase transition: Domain- 
wall renormalization-group studiesSi^ predict that the 
KGXY model in the strong disorder regime as well as the 
gauge-glass model undergoes a phase transition only at 
zero temperature. For the gauge-glass model, the zero- 
temperature transition was supported by computations 
of various quantities such as the current-voltage {IV) 
characteristics,^ the root-mean-square currentfi2*ii the 
correlation lengthjiSiH the glass susceptibilityjiSiiSii^ the 
autocorrelation functionfiiiH and the phase slip resis- 
tance^ A recent numerical consideration of the low- 
energy excitations also estimated Tc = in the gauge- 
glass modelii^ 

In contrast, the finite-size scaling analysis applied to 
the helicity modulus and the root-mean-square current 
in the RGXY model with strong disorder demonstrated 
that the superconducting phase transition of a non-BKT 
type occurs at finite Tj,, independent of the disorder 
strengthii In fact evidences for such a finite-temperature 
transition in the gauge-glass model have been presented 
in a few numerical studies computing the IV charac- 
teristics the correlation functionji the glass suscep- 
tibility,— and the linear resistance>ii Here it should be 
noted that critical exponents and Tc obtained from the 
resistively shunted junction (RSJ) simulations in large 
scaleai^ agree with other numerical studies of the gauge 
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glassjiii and lead to Tc = 0.22 (in units of Ej/ks, where 
Ej is the coupUng energy and ks the Bohzniann con- 
stant), the correlation length critical exponent u = 1.2, 
and the dynamic critical exponent z = 2.0. Furthermore, 
the barrier energy and the associated vortex mobility in 
the gauge glassi^ implies that superconducting order per- 
sists at low but finite temperatures. 

Very recently, an experiment has been performed on 
the PD JJAfiS and the critical temperature has been mea- 
sured as a function of the disorder strength, which re- 
veals a finite-temperature transition at strong disorder. 
Further, in the experimenfeiSi the scaling behavior of the 
IV characteristics is observed consistent with the corre- 
sponding numerical results for the gauge glassji^ii^ in- 
dicating the presence of a non-BKT type transition at 
a finite temperature in the strong disorder regime. A 
truly remarkable feature of this PDJJA experiment is 
that a very clear signal of a finite-temperature transition 
has been obtained for a single disorder realization of a 
200 X 800 square array. In contrast, typical attempts to 
settle the issue of the existence of a finite-temperature 
transition by direct simulations usually presume that a 
well converged disorder average is necessary. This in turn 
limits the simulations to small system sizes, so that any 
conclusion crucially depends on a priori assumptions for 
the finite-size scaling. The PDJJA experiment in Ref. 19 
suggests that one could avoid such a priori assumptions 
by instead obtaining well converged data for large sam- 
ples at the price of focussing on a single disorder realiza- 
tion. 

In the present paper we test this alternative strategy 
by performing RSJ simulations and computing the IV- 
characteristics of a large square array of Josephson junc- 
tions with one randomly chosen positional disorder re- 
alization. In parallel with the experimental results in 
Ref. Il9l we find strong numerical evidence for a finite- 
temperature transition. In view of the on-going contro- 
versy as to the existence of such a phase transition, we 
discuss the possible implications of our results. 

The present paper is organized as follows: In Sec. m 
we briefly describe the numerical method used in this 
work, and discuss the relation between the PDJJA and 
the KGXY model. Section llTTl presents mostly the results 
for the PDJJA at weak disorder, which is followed by 
the results at strong disorder in Sec. IIVI In Sec. we 
summarize the results of the present paper in relation 
with the existing studies. 



where (pij is the phase difference between the supercon- 
ducting islands at sites i and j. When an external mag- 
netic field B = Bz is applied, the magnetic bond angle 
Aij obtains the form 
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where a is the lattice constant and <I>o is the magnetic fiux 
quantum. Due to the positional disorder, the position of 
the zth island is given by 



{xi.yi) = {x°+5x„y^^+5yi) , 



(3) 



where = [x^, represents the ideal position without 
disorder and 5xi and 5yi are random quenched variables 
uniformly distributed in the interval [—A, A]. Hence- 
forth, Ti, Xi, 5xi^ and 5yi are all taken to be dimen- 
sionless, measured in units of the lattice spacing a. 

The magnetic frustration / is usually defined as the 
number of fiux quanta per plaquette. There has been 
intensive research interest in the role of the magnetic 
frustration in both classical^ and quantum^ systems. 
In the absence of positional disorder, every plaquette has 
an equal area, thus / is constant over the whole system. 
It is well known that the Hamiltonian |^ in this case is 
invariant under the transformation / / ± li2& In the 
present work, however, the plaquette area changes from 
place to place, and the Hamiltonian loses the above sym- 
metry, invalidating the equivalence between the cases / 
and /±1. We in this work define / as the disorder average 
(denoted to be [■■■]) of the flux through one plaquette 
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The magnetic bond angle Aij in Eq. I^Jl is correlated 
with the nearest neighboring bond angles since a change 
of position gives rise to changes of four magnetic bond 
angles A^ with site j being neighbors of site i. With- 
out such short-range correlations, Aij becomes a random 
quenched variable characterized by the disorder average 
[A,j\ = 27r/ (a;° -I- a;°) (?/° - y°). In case that / is an in- 
teger, one can gauge away the magnetic bond angle, es- 
tablishing the equivalence with the JiGXY model where 
the average frustration across the whole system vanishes. 
The variance of the sum of the magnetic bond angles 
around one plaquette in a PDJJA is given by 



II. 



MODEL AND SIMULATIONS 



With the assumption that the coupling energy is not 
affected by positional disorder, the Hamiltonian of the 
PDJJA in the presence of a transverse magnetic fleld 
readsS^ 



H = -Ej ^ cos ( 
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whereas the corresponding quantity in the HGXY model 
reads 
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where Aij has been taken to be uniformly distributed in 
[— TTT, m] with the disorder strength r. It is to be noted 
that unless A is close to unity, one obtains 

r « /A. (7) 

Consequently, the disorder strength of the PDJJA in 
comparison with the RGXY model is measured by /A, 
as found in previous studies of the PDJJA. ^^■^'^ It should 
be noted that the equivalence between the PDJJA and 
the KGXY model becomes valid only when the exter- 
nal magnetic field yields an integer value of the average 
magnetic frustration. The relation in Eq. ((TJ is practi- 
cally very useful since it provides a convenient way to 
perform experiments on the 2D KGXY model, in which 
the disorder strength can be tuned for one sample only 
by increasing the integral value of /. 

We below sketch briefly the numerical method adopted 
in this work. Introducing the twist variable D = 
(Dx,Dy) for the fluctuating twist boundary conditions 
(FTBC)^'' in L X i arrays, we write 0.^ in Eq. (Q) as 

0,, =e,-9,-rl--D, (8) 

where r°^- = — r° and 9i is the phase angle. Since the 
effects of positional disorder are introduced only through 
magnetic bond angles, we use the RSJ dynamics com- 
bined with the FTBC as follows: Using the local current 
conservation, the equation for the phase angle at site i is 
given by 

where time has been measured in units of h/2eicR with 
the single-junction critical current ic — 2eEj/h and 
the shunt resistance R (see Ref. for details), Gij is 
the square-lattice Green function, X]fc' denotes the sum- 
mation over the four nearest neighbor sites of j , and 
rjjk is the dimensionless thermal noise current satisfying 
iVi]) = and 

imj it)m (0)) - 2T [S.kSji - SaSjk) 5{t), (10) 

with the temperature T in units of Ej/ks- In order to 
use the efficient fast Fourier transform, we modify the 
periodic boundary condition for 6i according to 

Oi+Lii = + 2TifL5yi 

Qi+Ly = (11) 

where 2t: fLSyi is fixed in time and thus Oi-i^Ni = Oi- 
Here the usual periodic boundary conditions {Oi^^j^ = 
9j+Ly = Oi) are not applicable since Aij ^ Ai^^j^ j+L± 
for the PDJJA. 

We consider the system under uniform external cur- 
rents injected along the x-direction. The Josephson re- 
lation 2eVxlfi = 4>i+Nx i, where Vx is the voltage drop 



across the sample, together with the global current con- 
servation conditionpA leads to the equations of motion 
for the twist variables and Dy in the form 

Dx = sin {4>ij - Ajj) + riD^ - ig, 

= T2 H ^^'^ ^ ^^j ) + ' 

where y^/- \ denotes the directed sum over nearest 

neighboring bonds in the x{y) direction, is the external 
(driving) current density in units of the single-junction 
critical current ic, and rjo^ (jid ) is the thermal noise 
current for Dx {Dy) satisfying (rjo^) = {vOy) = and 

iVD. {t)VD. m = iVDy it)vDy (0)) = 2{T/L^)Sit). (13) 

We numerically integrate the equations of motion given 
by Eqs. ^ and ((T^ . and compute the IV characteristics 
with the average (dc) voltage V = Lv = —L{Dx), where 
V denotes the voltage drop per junction in units of icR- 

III. RESULTS FOR WEAK DISORDER 

Our simulation scheme requires sufficiently precise 
data for a large sample and a large disorder strength. We 
find that RSJ simulations of the PDJJA on a 128 x 128 
square lattice can be made to meet all these three re- 
quirements simultaneously, and thus use this lattice size 
in the present investigation. 

To probe phase transitions in weak and strong disorder 
regimes, we compute the IV characteristics at various 
temperatures and values of / with the latter controlling 
the disorder strength (see Sec.^. In our simulations we 
set the parameter for positional disorder equal to A = 
0.2. Since r « /A and the critical disorder strength in 
the RGXY model is believed to be Tc « OA^i^ we vary 
the average frustration / from 1 to 4 to cover both weakly 
(/ = 1 with /A « 0.2 < Tc) and strongly (/ > 2 with 
/A > Tc) disordered cases. 

Figure ^ presents the dc resistance v /id (in units of 
the shunt resistance R) versus the temperature T at var- 
ious frustrations / = 1,2,3, and 4. When / = 1, corre- 
sponding to the weakly disordered case, the resistance- 
temperature (RT) curves in Fig. disclose that the sys- 
tem undergoes a phase transition around T = 0.6, which 
is lower than Tc = 0.89 in a regular JJA.^** For strong 
disorder (/ > 2), phase transitions are observed to occur 
near T 0.2, hardly depending on the frustration. 

In order to determine the transition temperature in a 
more accurate way and to understand the dynamic criti- 
cal behavior in detail, we next investigate the IV charac- 
teristics at various values of / with the help of the scaling 
form suggested in Ref. |2^ In two dimensions, the cur- 
rent density and the electric field behave as J ^ T/^ and 
E ^ respectively, with the correlation length ^. 



FIG. 1; (Color online) Resistance vjid (in units of i?) versus 
temperature T (in units of EjjkB) for various frustrations 
in positionally disordered Josephson-j unction arrays. As the 
frustration is raised from / = 1 to / = 2, the effective disorder 
strength becomes larger, reducing the critical temperature. 
Larger values of / (= 3, 4) yield only insignificant variations 
of the critical temperature. 
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FIG. 2: (Color online) IV curves (v versus id) for / = 1 
in log scales at various temperatures T (in units of Ej Iks)- 
At r < 0.6, curves fit well to the power-law form v « i^^^, 
with the solid line representing the case z = 2. At higher 
temperatures, IV curves bend upwards, indicating a BKT- 
type transition. 

In the case Tc ^ 0, one replaces T in J ^ T/^ by Tc 
near criticality, and obtains according ly J - In this 
work, id and v correspond to the current density and the 
electric field, respectively, and yield the scaling form 

V - uC^F± (zrfO (14) 

with the scaling function F± above/below T^. At T — Tc, 
the correlation length ^ diverges and Eq. H14|) leads to 

In Fig. [3 we display the IV characteristics for / = 1, 
corresponding to the weak disorder regime, near T — 
Tc « 0.6 (see the RT curve for / = 1 in Fig.P). Due to 
the finite-size effects, IV curves exhibit the Ohmic behav- 
ior in the small current region as — > Oi^ The curves 
at r < 0.6 fit well to the power-law form w ~ z^: At 
T « 0.6, we have a ~ 3, and the lower T, the larger a. In 
contrast, the IV curves at higher temperatures (T > 0.6) 
clearly exhibit upward curvature, indicating the existence 



FIG. 3: (Color online) Slope d\nv/d\nid of the IV curve 
versus the driving current id at various temperatures (in units 
of Ej/kB). Due to finite-size effects, the slope first increases 
with id, then decreases toward the Ohmic value, unity, thus 
forming a peak. In a broad range of temperatures below Tc « 
0.7, the peak position i*i does not change substantially while 
it shifts to higher currents beyond Tc- Such behavior of the 
peak position is consistent with the BKT-type transition. 



of a finite current scale (and thus of a finite length scale) . 
From the relation a ~ z + 1 between the dynamic critical 
exponent z and the nonlinear IV exponent a, we thus 
reach the conclusion that a phase transition of the BKT 
nature occurs at Tc « 0.6 with z « 2. It should be noted 
that the scaling form in Eq. (|14|l is valid when the cor- 
relation length diverges only at Tc and becomes smaller 
as the temperature is varied (either lowered or raised) 
from Tc- Accordingly, when the phase transition is of the 
BKT type, characterized by the diverging length scale in 
the whole low-temperature phase, only high-temperature 
IV curves collapse to the scaling form, while the low- 
temperature part cannot be made to collapse due to the 
lack of the length scale ^. 

Alternatively, one can also plot the slopes of IV curves, 
given by dhiv/dlnid, versus the driving current id, as 
shown in Fig. |21 to confirm the BKT nature of the tran- 
sitioniSi As id is reduced, all curves should eventually 
crossover to the Ohmic behavior characterized by the 
unit slope {dlnv/dhiid = 1) due to the finite-size effects. 
The peak position of each curve in Fig. |3| measures 
a characteristic current scale, which is inversely propor- 
tional to the length scale in the system, i.e., ^ i/£,. 
The correlation length ^ increases as T approaches Tc 
from high temperatures. Near Tc in the high-temperature 
phase and also in the whole low-temperature phase, the 
correlation length becomes larger than the size of the 
system. In this case, the relevant length scale of the sys- 
tem is not the correlation length but the linear size L of 
the system, leading to ?2 ~ 1/L independent of the tem- 
perature. In Fig. 131 the peak position has almost the 
same value around 0.4 as T is increased from below up to 
T « 0.7; then appears to drift away toward larger val- 
ues as T is increased further beyond T — 0.7, which is in 
good agreement with what we expect for the BKT-type 
phase transition. We emphasize again that for a finite 
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system starts to move not at the true critical tem- 
perature but at a temperature higher than when the 
correlation length becomes comparable to L. Between 
T = 0.63 and 0.6, dlnw/rflni^ « 3 at ijj, which, com- 
bined with the power-law decay form in Fig. |3 makes 
us conclude that the BKT transition occurs at Tc ~ 0.6 
(in accord with the RT curve in Fig. ^ with the dy- 
namic critical exponent z w 2. This analysis of the IV 
curves clearly reveals that the phase transition of the 
BKT type occurs in the weak positional disorder regime 
of the PDJ JA, consistent with experimentiS. and previous 
numerical studies of JiGXY modeliS*^ 



IV. RESULTS FOR STRONG DISORDER 

In the strong disorder regime (/ > 2), we show in Fig.0| 
the IV curves for / = (a) 2, (b) 3, and (c) 4. It is clearly 
observed that for each value of /, there exists a well- 
defined temperature at which the IV curve shows a 
power-law behavior, manifesting the absence of a length 
scale at criticality. The critical temperature estimated 
from Fig. His = (a) 0.25, (b) 0.23, and (c) 0.23 for 
/ = 2,3, and 4, respectively, in agreement with Fig. ^ 
It is also shown in Fig. 0] that the IV curves at Tc fit 
well to the form i^^^ with the dynamic critical exponent 
z = 2. The IV curves in Figs. H (a), (b), and (c) aU 
bend upward above T^^ indicating that the PDJJA is in 
the high-temperature normal phase, while the opposite 
downward curvature below Tc implies the superconduct- 
ing phase in the limit of id 0. 

To find nature of the phase transition and the criti- 
cal temperature together with critical exponents in the 
strong disorder regime of the PDJJA, we employ the scal- 
ing form in Eq. (|14|l . which, with the correlation length 
^ - IT-Tc]-"", reads 

— — — ^ F± (iJT ~ Tc\-'') . (15) 

idlT-Tcl""" ^ ' ' ^ ' ' 

Figures El (d)-(f) show that data points in the IV char- 
acteristics [in Figs. 0] (a)-(c)] collapse into two function 
F± in Eq. 115|) : this confirms that for given value of 
/, there exists a finite-temperature phase transition of 
a non-BKT type. The critical temperature changes from 
Tc « 0.265 for / = 2 to Tc « 0.22 for / = 3 and 4, which 
is consistent with the experiment on the PDJJA^^ and 
numerical studies of the RGXY model"* and the gauge- 
glass modeliiS*iSiiLi^ On the other hand, in all three 
cases [/ = 2,3, and 4 in Figs. ^ (d)-(f)], the scahng 
collapse is achieved with v = 1.2 and z — 2.0, imply- 
ing that the nature of the transition in the strong disor- 
der regime remains unchanged as the disorder strength 
is increased. Furthermore, the dynamic critical exponent 
z — 2.0 agrees with the studies of dynamic behavior in 
the gauge-glass model^iiLii as well as with the experi- 
mental resultsi^. The critical exponent i/ = 1.2 obtained 
in this work is again consistent with the previous numer- 
ical resulta^iiSiii but not with the experimental result 



for the PDJJA, -'^^ the origin of which is not clear at this 
stage. 

We next investigate slopes of the IV curves for / = 4 
in the same manner as in Sec. Illll for the weak disorder 
case for / = 1. In Fig.jS] the slope d\nv/dlnid is plot- 
ted as a function of id at various temperatures around 
Tc- Note that finite-size effects yield the Ohmic behav- 
ior, dlnv/dhiid ~ 1, as id 0. It is observed that 
the curves at T = 0.23 and 0.20 are almost horizontal in 
broad ranges of the external currents, indicating scale- 
free behavior and in turn the existence of a phase transi- 
tion with the dynamic critical exponent z « 2.0, consis- 
tent with the finding in Fig.^f) for / = 4. Furthermore, 
the fact that the peak position changes rather abruptly 
from large currents to smaller currents at T = 0.23 to 
0.20 implies that the phase transition is of a non-BKT 
type, namely, the IV data are well scaled to two different 
functions, F± in Eq. (|15|l above/below Tc- 

In view of the ongoing controversy about the existence 
of a finite-temperature transition of a non-BKT charac- 
ter for strong disorder, one might ask what our results 
really imply. First of all, our evidence has been obtained 
for a finite sample with one disorder realization (albeit 
a very large one in comparison with those in most of 
earlier investigations). Does the phase transition survive 
in the large system size limit? Since a single positional 
disorder for a finite system can be periodically repeated, 
we can obtain an arbitrarily large system by just adding 
new squares with the same single disorder realization. 
Such an infinite system constructed from a single disor- 
der realization will, to our belief, most certainly have a 
finite-temperature transition of a non-BKT type. This 
means that the phase transition does exist per se. Sup- 
pose that we instead choose an arbitrarily large system 
and generate the disorder in the same random way as we 
have done for our 128 x 128 sample. Would the phase 
transition still survive? Here we find that various ran- 
domly generated disorder realizations for the 128 x 128 
system yield numerically very similar results, indicat- 
ing that for the size L = 128 there is already a large 
amount of disorder self-averaging. This again suggests 
that our results will survive in the large L limit, i.e., for 
the PDJJA with uniformly distributed disorder, chosen 
randomly, of strength A = 0.2. In order to check the 
finite-size effect on the self-averaging property in a more 
careful way, we have also computed voltages at T = 0.16 
in the strong disorder regime (/ = 4) for smaller sizes 
(L = 16, 32, 64, and 128) at two different disorder real- 
izations : The difference between voltages obtained from 
different disorder realizations is found to decrease as L 
is increased, and becomes negligibly small at L — 128. 
This indicates that indeed the self-averaging effect be- 
comes clear beyond L — 128. The next question is then 
whether our results also carry over to the KGXY model 
and the gauge-glass model. Here the evidence is more cir- 
cumstantial and based, on the one hand, on the strong 
connection between the PDJJA and the RGXY and ran- 
dom gauge-glass modelsiS, and, on the other hand, on the 
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FIG. 4: (Color online) IV characteristics {v versus id) at various temperatures for / = (a) 2, (b) 3, and (c) 4, and the 
corresponding scaling plots in (d)-(f). The solid lines in (a), (b), and (c) represent the power-law decay form v « i^^^ with 
z = 2. A well-defined temperature Tc separates the IV curves into two groups, one bending upwards and the other downwards. 
All data points in the corresponding IV characteristics are made to collapse into scaling functions, as shown in (d)-(f), with 
Tc = 0.265, 0.22, and 0.22, respectively. In all the three cases (d)-(f), the same critical exponents v = 1.2 and z — 2.0 are used, 
implying that the transitions in the strong disorder regime (/ = 2, 3, and 4) belong to the same universality class. 



strong similarity of the present results to some earlier re- 
sults obtained from the latter models i^iiLii 



V. SUMMARY 

The existence of a finite-temperature phase transition 
in the strong disorder regime of the PDJJA and the 
2D KGXY model, including the fully disordered case 
of the 2D gauge-glass model, are still in an intensive 
debate. This work has been motivated by the very re- 
cent experiment on 800 x 200 Josephson junction arrays 
with positional disorder, and explored numerically the 
dynamic critical behavior of the PDJJA in transverse 



magnetic fields. Adopting the RSJ dynamics,^"* we have 
computed the IV characteristics of the PDJJA with the 
positional disorder parameter A = 0.2 for frustration 
/ = 1,2,3 and 4. The relation r w /A between the 
disorder strength r in the 2D KGXY model and the pa- 
rameter A in the PDJJA, combined with the critical dis- 
order strength Tj. « 0.4)2ii4 implies that the PDJJA for 
f ~ 1 corresponds to the weakly disordered case while 
strong frustration (/ > 2) puts the PDJJA in the strong 
disorder regime. The scaling analysis^® of IV curves and 
their slopes^^ has revealed clear evidence for Tc 7^ in 
the PDJJA with strong disorder, which agrees with the 
experiment on the PDJJA^^ as well as previous numer- 
ical studies of the KGXY model^ and the gauge-glass 



7 



4.5 
4 
3.5 
f 3 
3 2.5 
2 
1.5 
1 



r = 0.16 
r = 0.18 

T = 0.20 
T = 0.23 ^— 
r = 0.25 
r = 0.28 

r = 0.30 --B-- 

T = 0.40 
T = 0.50 

r = 0.60 --+-- 




model 



7.15.16.17.18 



0.01 



0.1 



FIG. 5: (Color online) Slope {d\nv/d\nid) of IV curves 
against id for f — 4, corresponding the strong disorder regime 
of the PDJJA, at various temperatures T (in units of Ej/kB)- 
Power-law behavior appears at T = 0.23 to 0.20 with the dy- 
namic critical exponent z ~ 2, indicating the existence of a 
phase transition with a diverging length scale only at Tc- 
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Figure exhibiting the phase diagram on the plane 
of the temperature and frustration summarizes the re- 
sults of the present work. The critical temperature Tc is 
observed to reduce rapidly as the average frustration /, 
controlling the disorder strength, is increased in the weak 
disorder region. It then appears to saturate toward a con- 
stant value Tc = 0.22 in the strong disorder regime. Be- 
low the phase boundary separating the superconducting 
phase at low temperatures and the normal phase at high 
temperature, there exist two different superconducting 
orders, according to the disorder strength: In the weak 
disorder regime, e.g., f ~ 1, the superconducting state 
is the low-temperature BKT phase characterized by the 
divergence of the correlation length. On the other hand, 
for strong disorder (/ ^ 2), the transition is of the non- 
BKT type with the well-defined correlation length criti- 
cal exponent i' = 1.2, consistent with the value obtained 
previously for the 2D gauge-glass modeliiiSiii but incon- 
sistent with the experimental finding i' = 2.0±0.3.^'^ The 
origin of this discrepancy is not clear at present and needs 
more detailed investigation in the future study. Finally, 
we point out that the resemblance between the phase dia- 
gram in Fig. inland the corresponding diagram obtained in 
the recent experiment on the PDJJA^^ is striking, which 
is also in very good agreement with the numerical study 
of the 2D random gauge XY modeli^ 



FIG. 6; Phase diagram on the plane of temperature T (in 
units of Ej/kn) and frustration / (controlling the disorder 
strength). The solid line, separating the high-temperature 
normal phase from the low-temperature superconducting 
phase, is only a guide to the eye. The critical temperature 
saturates to a finite nonzero value Tc ~ 0.22 as / is increased. 
Below Tc, depending on the disorder strength, there exist two 
different superconducting phase: a BKT-like phase for / = 1 
and a non-BKT type one for / > 2. 
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